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Two-dimensional (2D) semiconducting transition metal dichalcogenides (TMDs) can withstand a 
large deformation without fracture or inelastic relaxation, making them attractive for application 
in novel strain-engineered and flexible electronic and optoelectronic devices. In this study, we 
characterize the mechanical response of monolayer group VI TMDs to large elastic deformation 
using first-principles density functional theory calculations. We find that the ultimate strength and 
the overall stress response of these 2D materials is strongly influenced by their chemical composition 
and loading direction. We demonstrate that differences in the observed mechanical behavior can 
be attributed to the spatial redistribution of the occupied hybridized electronic states in the region 
between the transition metal atom and the chalcogens. In spite of the strong covalent bonding 
between the transition metal and the chalcogens, we find that a simple linear relationship can be 
established to describe the dependence of the mechanical strength on the charge transfer from the 
transition metal atom to the chalcogens. 



I. INTRODUCTION 

Since its first mechanical exfoliation from graphite, 
graphene — an allotrope of carbon with a two- 
dimensional (2D) structure — has attracted tremendous 
scientific and technological attention due to novel elec- 
tronic, mechanical, and chemical properties!^ The in- 
tense interest in graphene has also stimulated an active 
search for possible inorganic 2D or quasi-2D materials 
with unique characteristics. Examples of such mate- 
rials include hexagonal boron nitride, transition metal 
dichalcogenides (TMDs), and transition metal oxidesP^ 
Earlier attempts to obtain monolayer and a few-layer 
thick samples of these materials were based on mechani- 
cal exfoliation from their bulk counterparts. Recent ad- 
vances in fabrication techniques, such as liquid exfolia- 
tion, lithium intercalation, epitaxy, laser thinning, and 
chemical vapor deposition, have successfully resulted in 
a scalable production of large 2D samples containing only 
a single layer to a few layers PUGU 

Among the non-carbon 2D materials that are being 
explored, dichalcogenides of group IV, V and VI transi- 
tion metals are of particular interest. Similar to graphene 
and its bulk counterpart graphite, these materials have 
strong in-plane bonding while the individual layers are 
bonded by weak van der Waals interactions. Several 
TMDs, unlike graphene, have an intrinsic band gap in the 
range of 1 - 2 eV, making them attractive for field effect 
transistors and optoelectronic devices 1121121 Many single- 
layer TMDs also demonstrate novel electronic properties 
that are otherwise not seen in their bulk or few-layer 
thick counterparts, for example, direct to indirect band 
gap transition, enhanced photoluminescence, and valley 
p ol ar iz at ion f^EU 

Recent mechanical experiments show that 2D TMDs 
can sustain very large elastic strains (~ 10% effective 



in-plane strain) with a high resistance to inelastic relax- 
ation and fracture, and hence these materials are also of 
interest for flexible electronics applications. 22 Since the 
electronic conduction and valence energy states both de- 
pend on the strain, it has also been proposed that elastic 
strain can be employed to tune the band gap in mono- 
layer TMDsP^ID In order to realize the potential ap- 
plications of various 2D TMDs in strain-engineered and 
flexible devices, an accurate characterization of their abil- 
ity to withstand large elastic deformation is crucial. 

Recently, Bertolazzi et alW^ performed the experimen- 
tal measurements on the stiffness and breaking strength 
of monolayer M0S2 by using an atomic force micro- 
scope tip to deform the monolayer M0S2 placed on a 
prepatterned Si02 substrate. They found that the ef- 
fective Young's modulus and average breaking strength 
of monolayer MoS 2 are 270 ± 100 GPa and 23 GPa, re- 
spectively. Using a similar technique, Castellanos-Gomez 
et a/P^ studied the mechanical properties of freely sus- 
pended M0S2 nanosheets with 5 to 25 layers and ob- 
tained the mean Young's modulus of 330 ± 70 GPa. 
These experiments clearly establish 2D M0S2 as an ultra- 
strength material able to withstand large elastic defor- 
mation. However, such studies provide only effective 
orientation-averaged values of the mechanical properties 
and are limited only to single composition, while the rest 
of the 2D semiconducting TMDs remain yet to be charac- 
terized. The theoretical studies of mechanical properties 
of monolayer TMDs are also limited to the calculation of 
simple elastic properties that are relevant only at small 
deformations) 2 ^ 2 ^ What is the intrinsic ultimate strength 
— the maximum stress a defect-free material can with- 
stand without failure — of monolayer TMDs along dif- 
ferent crystal directions and how is it dependent on the 
chemical composition? Indeed, the relationship between 
the intrinsic stress response of monolayer TMDs to large 
elastic deformation, the atomic-level geometry, and the 
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chemical composition remains unclear. 

In this article, using first-principles density functional 
theory simulations, we report on the fundamental me- 
chanical behavior of monolayer TMDs when subjected to 
a large deformation. Here we consider TMDs with chem- 
ical composition MX2, where M (= Mo, W) is a group 
VI transition metal and X (= S, Se, Te) is a chalcogen. 
Among the range of naturally occurring TMDs, group VI 
TMDS are semiconductors and, therefore, are good can- 
didates for flexible and strain-engineered optoelectronic 
and photonic devices. Here we have determined the en- 
tire stress response of monolayer TMDs to large applied 
elastic strains — the intrinsic materials characteristics 
that are crucial in designing the strain-engineered and 
flexible devices. We find that the stress response and 
the ultimate strength of monolayer TMDs strongly de- 
pends on the chemical composition as well as the load- 
ing direction. For the same transition metal M, the sul- 
fides (MS2) are the strongest and selenides (MSe2) are 
stronger than the tellurides (MTe2). The chalcogens of 
W (WX2) can accommodate larger stresses than that of 
Mo (M0X2). For all compositions, the ultimate strength 
of the monolayer along the armchair direction far exceeds 
the strength in zigzag direction, typically by a factor of 
1.5 - 2.2. Our electronic structure analysis reveals that 
the origin of the observed mechanical behavior can be 
attributed to the hybridization between M-d and X-p or- 
bitals and the charge transfer from the transition metal 
to the chalcogens. 



II. RESULTS AND DISCUSSION 

All MX2 TMDs monolayers considered in our study 
share a similar hexagonal crystal structure as shown in 
Fig. [l] A single layer of MX2 consists of three X-M- 
X sublayers, with hexagonal lattice plane of the transi- 
tion metal M sandwiched between identical hexagonal 
planes of the chalcogens X. In order to evaluate the 
mechanical response of the pristine TMD monolayers, 
we carried out first-principles density functional calcu- 
lations within the local density approximation as im- 
plemented in ABINiTpH] 

Norm-conserving pseudopo- 
tentials in the form of Hartwigsen-Geodecker-Hutter 
pseudopotentiald^l including the semicore electrons as 
valence electrons for Mo and W were used to describe 
the interaction between core and valence electrons. An 
energy cutoff of 65 Hartree was used for plane wave basis 
expansion. The fc-point sampling schemes of 15 x 15 x 1 
and 10 x 15 x 1 were employed during the geometry op- 
timization of strain-free monolayer sheets and the calcu- 
lation of stress-strain relations, respectively. 

For each of the MX2 monolayers, we first performed 
the structure relaxation including the lattice constants 
and atomic coordinates as summarized in Table |U The 
lattice constants obtained from our DFT calculations are 
in excellent agreement with the experimental measure- 
ments (error less than 1%). 12 Our results are also in close 




FIG. 1. A schematic showing the crystal structure of a TMD 
monolayer (MX2) with (a) top view and (b) side view. The 
transition metal (M) and chalcogen (X) atoms are represented 
by purple (large) and yellow (small) spheres, respectively. The 
unit cell (red) and the orthogonal supercell (blue) are also 
depicted in (a). The armchair and zigzag orientations corre- 
spond to x and y directions, respectively. 



TABLE I. Lattice constants, M-X bond lengths and X-X dis- 
tances for monolayer MX 2 (M = Mo, W; X = Se, Se, Te) 
TMDs. 



MX 2 


a (A) 


M-X (A) 


x-x (A) 


MoS 2 


3.13 


2.38 


3.11 


MoSe 2 


3.25 


2.51 


3.32 


MoTe 2 


3.49 


2.70 


3.60 


WS 2 


3.13 


2.39 


3.12 


WSe 2 


3.25 


2.51 


3.33 


WTe 2 


3.49 


2.70 


3.61 



agreement with those computed with different methods 
such as numerical atomic orbit alsp^ PAW potentials^. 

Similar to the 2D hexagonal lattice of graphene, two 
high-symmetry directions can be identified in the crys- 
tal structure of monolayer TMDs, namely, armchair and 
zigzag directions, which are oriented along and perpen- 
dicular to M-X bonds when projected on the plane con- 
taining transition metal atoms, respectively. Since the 
mechanical behavior of a defect-free material is ulti- 
mately controlled by the strength of its chemical bonds, it 
can be expected that the stress response of the monolay- 
ers along the armchair and zigzag directions will also be 
distinct. We therefore employed an orthogonal supercell 
as shown in Fig. [IJa) to calculate the stress-strain curves 
along the two high symmetry directions. In each case, the 
super-cell was deformed along the high symmetry direc- 
tion in small strain increments and the in-plane supercell 
vector normal to the the applied strain was allowed to 
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FIG. 2. Tensile stress a as a function of uniaxial strain e along (a) armchair and (b) zigzag directions, respectively, for monolayer 
MX 2 (M = Mo, W; X = S, Se, Te) TMDs. Solid and dashed lines are used for WX 2 and MoX 2 , respectively. 



relax in order to account for the Poisson's contraction 
as typically observed in the experiments. The super cell 
stress as output directly by ABINIT is the stress aver- 
aged over the supercell volume. We renormalized the 
calculated supercell stress by a factor of Z/h, where Z is 
the vacuum distance 12 A used to avoid the interaction 
between periodic images and h is the interlayer distance 
in the bulk MX 2 . The strain is denned as e = 
where Lq and L are the supercell length along armchair or 
zigzag directions before and after applying tensile strain, 
respectively. 

In Figs.[2|a) and[2|b) we present the calculated stress- 
strain relations along armchair (x) and zigzag (y) direc- 
tions, respectively. It can be readily noted that for small 
strain, the stress for all MX2 exhibits linear dependence 
on the applied strain for both loading directions. Ta- 
ble |TI| presents the Young's modulus E calculated as the 
slope of the stress-strain curve in the region of strain less 
than 4%. It is evident that for all compositions, the val- 
ues of the Young's modulus obtained independently from 
the stress-strain curve for armchair and zigzag loading 
direction are virtually identical. This behavior is consis- 
tent with the symmetry of the MX2 crystal lattice - the 
Young's modulus and other second order elastic constants 
are essentially isotropic due to the hexagonal symmetry 
in the basal plane. Our calculated values for the Young's 
modulus are also in good agreement with earlier studies. 
For instance, we obtained the Young's modulus of M0S2 
to be ~ 220 GPa, while experimental measurements have 
reported values of 270 ± 100 GPa and 330 ± 70 GPa for 
monolayer ^ and few layer M0S2 sheets^, respectively. 

As MX2 monolayers are strained further (e > 4%), 
the stress-strain response deviates from the linear behav- 
ior. For large strains, the hexagonal symmetry is broken, 
with the stress developed upon loading in armchair di- 
rection (Fig. |2|a)) much larger than in the zigzag di- 
rection (Fig. f2Tb)) . Upon straining further, the stress 
continues to increase till it reaches a maximum, termed 



as the ultimate strength a*. Table [TT] presents the cal- 
culated values for the peak stress a* and the corre- 
sponding ultimate strain £* for all chemical composi- 
tions. It can be observed that in general, the chalcogens 
of W (WX2) have larger moduli and tensile strength than 
that of Mo (M0X2). For the same transition metal, sul- 
phides (MS2) are the strongest while tellurides (MTe2) 
are weakest. Moreover, Table [TT] also presents values 
for the anisotropy factor <p = Car/^zz? which measures 
the relative strength along the two high symmetry direc- 
tions. We find that the anisotropy in stress response is 
inversely correlated with the strength of the monolayer 
sheets — the MX2 with lower Young's modulus and ulti- 
mate strength (for example, tellurides) are characterized 
by larger anisotropy factors. 

Our calculated value for the average ultimate strength 
of monolayer M0S2 (22 GPa) is in excellent agreement 
with the value of 23 GPa obtained by the atomic force 
microscopy measurements. Furthermore, the observed 
stress-strain response of the single-layer structure of 
the TMDs can be compared with their nanotube coun- 
terparts. For instance, density functional tight bind- 
ing (DFTB) simulations have predicted a Young's modu- 
lus of 209.7 GPa and 236.6 GPa for armchair and zigzag 
nanotubes of comparable diameter about 12.0 A, respec- 
tively, in close agreement with our results. 33 Moreover, 
their ultimate strength was reported to be 29.1 GPa and 
31.6 GPa, respectively, again in good qualitative agree- 
ment with our results. However the DFTB studies of 
M0S2 nanotubes suggested similar value (~ 16%) for the 
ultimate strain £* corresponding to the ultimate strength 
for both armchair and zigzag directions, in contrast to 
our results where much bigger values of e* (> 28%) are 
observed for the armchair loading direction compared to 
e* (< 19%) for the zigzag loading direction (see Table [ill). 

It is also worth noting that a similar transition in 
stress-strain curve from an isotropic and linear stress re- 
sponse at small strains to anisotropic and nonlinear stress 
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TABLE II. Calculated Young's modulus E (GPa), ultimate strength cr* (GPa), ultimate strain £* corresponding to the ultimate 
strength for armchair and zigzag directions, anisotropy factor <j> and the charge transfer AQ (e) from M to X for strain- free 
MX2 monolayer sheets. 



MX 2 


E 


Armchair (x) 

* 

a 


e* 


E 


Zigzag (y) 
cr* 


e* 





AQ 


M0S2 


222.75 


27.35 


0.28 


219.46 


16.90 


0.19 


1.62 


0.92 


MoSe 2 


178.78 


22.68 


0.29 


175.97 


12.86 


0.16 


1.76 


0.73 


MoTe 2 


125.94 


17.12 


0.32 


123.54 


7.88 


0.14 


2.17 


0.37 


WS 2 


244.18 


29.96 


0.28 


240.99 


19.91 


0.18 


1.50 


1.07 


WSe 2 


196.81 


24.70 


0.30 


194.13 


15.05 


0.17 


1.64 


0.83 


WTe 2 


137.32 


18.71 


0.32 


135.27 


9.30 


0.15 


2.01 


0.44 
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FIG. 3. Total and projected density of states for strain-free monolayer MX 2 (M = 
transition metal M and the chalcogen X, the density of states are projected onto s, 
Fermi level is set to zero. 
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response at large strains has also been observed in the 
case of gr ap he neM^ However, there is one remarkable dif- 
ference: while TMDs and graphene both have hexagonal 
lattice, they exhibit a contrasting anisotropy in stress re- 
sponse. The peak stress in graphene for armchair tension 
is 9% lower than that for zigzag direction, in contrast to 
the behavior seen in Fig. [2] 

In order to identify the atomic-level origin of the me- 
chanical behavior of various transition metal dichalco- 



genides considered here, we analyzed their electronic 
structures by using the Vienna ah initio simulation pack- 
age (VASP) 35 with LDA exchange correlation functional. 
Frozen-core projector-augmented wave (PAW) method 
was employed to describe the interaction between core 
electrons and valence electrons.^ An energy cutoff of 400 
eV was used for plane wave basis expansion. 

The calculated total and projected density of states of 
unstrained MX2 monolayer sheets are presented in Fig. [3] 



Since the strength of the interatomic bonding and the 
consequent mechanical response of the materials to defor- 
mation are determined primarily by the occupied states 
just below Fermi level, we focus on the energy range - 
8 eV below the Fermi level. It can be seen that in this 
range, the outermost d orbitals of the transition metal M 
overlap significantly with the outermost p orbitals of the 
chalcogen X, indicating a strong p-d hybridization. Fur- 
thermore, a careful examination of the peaks of these hy- 
bridized states yields a crucial insight on the composition 
dependence of the mechanical response seen in Fig. [2] For 
instance, the peaks of p-d orbitals in M0S2 are located 
at -5.25, -4.02, -3.27, -1.93 and -0.83 eV below the Fermi 
level. In MoSe2, these peaks shift toward the Fermi level 
by 0.28, 0.13, 0.15, 0.03 and 0.11 eV, respectively. The 
peaks shift further toward the Fermi level in MoTe2 by 
0.31, 0.52, 0.78, 0.30 and 0.08 eV, respectively. In addi- 
tion to the change in peak positions, the height of peaks 
also changes as the composition is changed from M0S2 
to MoTe2i the peaks of low- lying Mo-4<i states are low- 
ered while the peaks of states close to Fermi level are 
enhanced. As seen in Fig. |3j the hybridized states for 
the chalcogens of W follow an essentially similar trend. 

The changes in the p-d hybridized orbitals in monolayer 
TMDs can be quantitatively characterized by computing 
the center of d bands below Fermi level, which is defined 



as 



E d 



fz F PDOS(£,d) x EdE 



fz F PDOS(£,d) dE 



where PDOS(£^, d) is the density of states projected onto 
the d orbitals of transition metals Mo and W, the lower 
end energy El is taken as -8 eV to cover the p-d hy- 
bridized states and the Fermi level Ep is set to zero. 
For M0S2, MoSe 2 and MoTe 2 , the calculated d band 
centers Ed are -2.68, -2.37, -1.94 eV, respectively. For 
WS 2 , WSe 2 and WTe 2 , E d is located at -2.98, -2.65, - 
2.18 eV, respectively. Since the hybridization of p-d or- 
bitals indicate a covalent electron sharing between M and 
X atoms, the deeper the center of the d bands below the 
Fermi level, the stronger is the M-X bond. Consequently, 
the tellurides with shallow E d demonstrate lower Young's 
modulus and tensile strength compared to sulphides and 
selenides. Chalcogens of W have deeper E d than the 
chalcogens of Mo and are therefore stronger. The analysis 
of the locations of d band centers presented here explains 
the observed composition dependence of the mechanical 
response of the monolayer TMDs as seen in Fig. |2j 

In order to visualize the hybridized electronic states be- 
tween the transition metals and chalcogens, we have fur- 
ther calculated the electronic charge distribution. Fig. [4] 
shows the bonding charge density in the plane passing 
through both the transition metal and chalcogen atoms. 
The bonding charge density is obtained as the difference 
between the valence charge density of strain-free MX2 
sheet and the superposition of the valence charge density 
of the constituent atoms. A positive value (red) indicates 
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FIG. 4. Bonding charge density for monolayer MX2 (M = 
Mo, W; X = S, Se, Te) TMDs, obtained as the charge density 
difference between the valence charge density of the monolayer 
and the superposition of the valence charge density of the 
neutral constituent atoms. Red and blue colors indicate the 
electron accumulation and depletion, respectively. The color 
scale is in the units of e/Bohr 3 . 
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FIG. 5. Variation of (a) the Young's modulus and (b) the 
ultimate strength of monolayer MX 2 (M = Mo, W; X = S, 
Se, Te) TMDs with the charge transfer AQ from transition 
metal M to chalcogens X. 
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TABLE III. Fitted values in the linear relation ~ a + bAQ 
to describe the dependence of the Young's modulus E and 
ultimate strength cr^ R and cr^z on the charge transfer AQ. 





a 


b 


E 


171.97 


61.34 


_* 


18.21 


10.40 


* 


16.86 


1.60 



electron accumulation while a negative value (blue) de- 
notes electron depletion. These bonding charge distri- 
butions clearly show the electron accumulation in the 
middle region of M and X. The amount of charge local- 
ized in this region qualitatively indicates the strength of 
the covalent M-X bond. It is evident that as the com- 
position changes from sulphides to tellurides, the charge 
accumulation in the bonding regions gradually become 
less, indicating a weakening of M-X bond, lower Young's 
modulus and ultimate strength. 

Fig.[4]also shows a large amount of charge transfer from 
the transition metal to the chalcogens. This charge redis- 
tribution can be quantitatively estimated by computing 
the charge transfer from M to X. Table [IT] presents the 
magnitude of t he ch arge transfer obtained using Bader 
charge analysis We depict the Young's modulus E 
and ultimate strength a* as a function of charge trans- 
fer AQ in Figs. 5(a) and |5(b)| respectively. It is evi- 
dent that the mechanical properties of transition metal 
dichalcogenides exhibit a linearly increasing relation with 
the charge transfer from transition metal atom to chalco- 
gen atoms. By fitting the values of the Young's modulus 
and the ultimate strength with the charge transfer, this 
relationship can be simply described as E(a*) ~ a + bAQ 
where the values of fitted parameters a and b are listed 



in Table III These expressions provide an approximate 
and simple description of the mechanical properties of 
transition metal dichalcogenide monolayers sharing sim- 
ilar crystal structures. 



group VI TMDs to large elastic deformation, the atomic- 
level structure, and chemical composition. Our calcu- 
lations demonstrate that the chemical composition of a 
monolayer TMD strongly influences its mechanical re- 
sponse to large elastic deformation, with the chalcogens 
of W exhibiting much greater ultimate strength than 
the chalcogens of Mo. The stress response of monolayer 
TMDs also depends on the crystal symmetry, with the 
armchair directions being consistently stronger across all 
chemical compositions. The origin of the observed me- 
chanical behavior can be attributed to a strong hybridiza- 
tion between the outermost p orbitals of the chalcogens 
and the d orbitals of the transition metal. This hybridiza- 
tion leads to a redistribution of the electronic charge 
to the shared region between the transition metal and 
chalcogen atoms. Our study clearly highlights the inter- 
play between the atomic structure and the composition 
of monolayer TMDs that lends them a great strength 
to sustain large reversible deformations — deformations 
large enough to modify their electronic structure as en- 
visioned for novel optoelectronic and photonics applica- 
tions. It should be emphasized that our first-principles 
investigation provides a realistic estimate of the intrinsic 
mechanical response of monolayer TMDs as measured by 
the nano-indentation experiments. In such experiments, 
the stressed region of the monolayer under the atomic 
force microscope tip can be expected to be defect-free 
and, therefore, the measured stresses can approach the 
ultimate strength. Therefore, our investigation of the 
stress-strain relation of monolayer TMDs presented here 
not only yields crucial insights into mechanical behav- 
ior of these materials but also can be compared with the 
experimental measurements. 
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